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ABSTRACT 



Recently, Park & Gott (1997) claimed that there is a statistically significant, strong, 
negative correlation between the image separation A9 and source redshift z s for grav- 
itational lenses. This is somewhat puzzling if one believes in a flat (fc = 0) universe, 
since in this case the typical image separation is expected to be independent of the 
source redshift, while one expects a negative cor relation in a k = — 1 universe and a 
positive one in a fc = +1 universe. Park fc Gott explored several effects which could 
cause the observed correlation, but no combination of these can explain the observa- 
tions with a realistic scenario. Here, I explore this test further in three ways. First, I 
show that in an inhomogeneous universe a negative correlation is expected regardless 
of the value of k. Second, I test whether the A9-z s relation can be used as a test to 
determine Ao and Qq, rather th an just the sign of k. Third, I compare the results of 
the test from the Park fc Gotl] sample to those using other samples of gravitational 
lenses, which can illuminate (unknown) selection effects and probe the usefulness of 
the A8-z s relation as a cosmological test. 

Key words: gravitational lensing - cosmology: theory - cosmology: observations. 



1 INTRODUCTION 

Historically, there has been little interest in the A9-z s rela- 
tion compared to other cosmological tests based on grav- 
itational lensing sta tistics, perh aps because the inflation- 
ary paradigm (e.g. 



Guth 1981), which began about the 



s ame time as the discovery of the first gravitational lens 



rod') relations, 'conventional' gravitational lensing statis- 
tics, age of the universe) are relatively insensitive to the ra- 
dius of curvature of the universe (Ro ~ (|fio + -V) — l|) _5 )i 
being degenerate in combinations of Ao and Qo in direc- 
tions roughly perpendicular to lines of constant Ro in the 
Ao-fio plane. A notable exc eption are constraints de r ived 



(Walsh) Carawoll fc Woy i 
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haa be 
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from CMB anisotropic s (e.g. fccott, White fc Silk 1995 



ential. Since a flat (fc — 0) universe is a robust prediction of 



Sugiyama fc Silk 1997) 



Hu 



inflation, many researchers assume this and consider only 
flat universes (or, at most, k = — 1 cosmological models 
with Ao = 0). Due to the fact that for the popular singular 
isothermal sphere model for a single-galaxy lens the aver- 
age image separation A9, integrated over the lens redshift 
Zd from 2d = to Zd = z s , is completely independent of the 
source redshift z s in a flat universe, there is little point in 
pursuing the A9-z s relation if one is interested primarily in 
flat cosmological models. If one is not committed to a flat 
universe, then of course one should not assume k = 0, but 
even if one believes that the universe must be flat, it is still 
important to test this belief observationally. The situation is 
somewhat worsened by the fact that most 'standard' cosmo- 
logical tests such as the m-z (magnitude-redshift or 'stan- 
dard candle') and 9-z (angular size-redshift or 'standard 
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2 THEORY 

For a singular isothermal sphere lens, the angular image sep- 
aration is given by (e.g. Turner, Ostriker & Gott 1984) 



—art • 



(i) 



where v is the velocity dispersion and D is the angular size 
distance (see below). Even if the singular isothermal sphere 
is not a perfect model for the gravitational lens systems 
considered, it is still a good approximation when one is con- 
cerned only with the image separation. For a given v , by 
combining Eqs. (5) and (6) in Gott, Park fc Lee (1989 ) and 



using the more appropriate and more general angular size 
distances, one obtains an expression for the average image 
separation A9, by integrating over the lens redshift Zd from 
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Zd = to Zd = z B 



A9{z B 



dz d 



DlDj (1 + Zdf 
DlQ 



A9(0) 



(2) 



dz d 



D 2 ds D 2 d (l + z d f 
DlQ 



where 



Q 



^n (l + z d f - (fio + A - 1) (1 + z d f + Ao 



(3) 



The Dij (with := Dot) in Eqs. and (g) are angular 
size distances, which are functions of the lens and source red- 
shifts Zd and 2 a , the cosmological parameters Ao and ilo as 
well as the 'homogeneity parameter' 77, which gives the frac- 
tion of smoothly, as opposed to dumpily, distributed matter 
along the line of sight. Note that Eq. (m) is valid for all com- 
binations of Ao, fio and 77. The angular size distances can be 
computed for arbitrary combinations of these parameters b y 
the method outlined in [Kayser, Helbig fc Schramm (1997| ). 

Figures [j] and ^| show A9 as a function of z B for var- 
ious cosmological models, for 77 — 1 (the traditional case 
assuming a completely homogeneous universe) and 77 = as 
extreme cases. Note in Fig. hi that the curve is a horizontal 
line for k — 0, has positive slope for k = +1 and negative 
slope for k = — 1, where k := sign(fio + Ao — 1). In Fig. H, for 
77 = 0, the slope is negative regardless of the value of k. Thus, 
at first sight it appears that an inhomogencous u niverse, a 
possibility not investigated by Park & Gott (1997, hereafter 
PG), might be able to explain the puzzling negative correla- 
tion between A9 and z B . However, it is shown in Sect. |E]that 
even the extreme 77 = scenario produces an anticorrelation 
which is much weaker than that found by PG. This effect 
can be qualitatively understood by realizing how Eq. (^) is 
affected by decreasing 77: inspection shows that this might 
be estimated by examining Dd B /D a . All other things being 
equal, the angular size distance increases with decreasing 77. 
Also, the effect of 77 is more noticeable at large redshift dif- 
ferences. Since z B > z B — Zd, the denominator is the more 
important term, and so decreasing 77 increases D a and so 
decreases Dd B /D B and thus A9(z B )/ A9(0). 



3 DATA 

PG used an inhomogeneous sample of gravitational lenses 
from the literature. While this seems problematic at first 
sight, PG noted that there is no reason to believe that this 
should influence the analysis. Nevertheless, it is worth com- 
paring the PG results to those obtained from a better defined 
sample. 

The observational data provided by the JVAS and 
CLASS surveys offer an independent sample of gravitational 
l enses. JVAS is the J odrell Bank VLA Astrometric Survey 
( IPatnaik et al. 1992ft; CL ASS is the Cosmic Lens All-Sky 
Survey ( |Mycrs ct al. 199*^ ) . Even though the observational 
tasks are not yet complete, the JVAS and CLASS surveys 
which constitute the database have already yielded sufficient 
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Figure 1. Normalized image separation as a function of source 
redshift. From the top, the (Ao,f2o) values are (2,4), (0,4), k = 0, 
(0,0.7), (0,0.3) and (-5,1). For k = the result is valid for all 
(Ao,f2o) values whose sum is 1. 77 = 1. 
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Figure 2. The same as Fig. hi except that here 77 = 0. 



gravitational lenses to enable one to make an independent 
analysis. Table [l] shows the current state of knowledge about 
the JVAS/CLASS gravitational lenses. Note that the ques- 
tionable source redshift for 2114 + 022 is probably the red- 
shift of an additional lensing galaxy (this interpretation is 
supported by several independent lines of evidence). 

Although not all source redshifts in the JVAS/CLASS 
sample are known, 8 out of 11 are, and based on our sur- 
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Table 1. The JVAS/CLASS gravitational lenses. 



Name 


# images 


A(9 

[arcsec] 


lens galaxy type 


2d 


2s 


0218+357 


ring + 2 


0.33 


spiral 


0.6847 


0.96 


0414+0534 


4 


2.0 


elliptical 


? 


2.62 


0712+472 


4 


1.2 


? 


0.406 


1.339 


1030+074 


2 


1.6 


peculiar 


0.599 


1.535 


1422+231 


4 


1.2 


? 


0.65 


3.62 


1600+434 


2 


1.4 


spiral 


0.4144 


1.589 


1608+656 


4 


2.2 


spiral? 


0.64 


1.39 


1933+503 


4+4+2 


0.9 


? 


0.755 


? 


1938+666 


4+2 


0.9 


? 


? 


? 


2045+265 


4+1? 


2.0 


? 


0.87 


1.28 


2114+022 


2+2? 


2.4 


? 


0.316 


0.588? 
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Figure 3. Source redshifts z s and image separations A8 (in 
arcsec) for the gravitational lens systems studied in this paper. 
Crosses represent the PG sample (20 systems; note that two data 
points with A9 £s 6 arcsec almost coincide); diamonds represent 
the JVAS/CLASS sample (8 systems; of course only those with 
known source redshifts are included). Note that there is an over- 
lap of four data points. The filled diamond represents the system 
0218 + 357, which was not used by PG although its sour ce redshift 
had b een published before the PG analysis was done (Lawrence 
1996|). 



vey, discovery and followup strategies there is no reason 
to suspect the unknown source redshifts to be statistically 
different from those already known. Figure ^| shows the 
source redshifts and image separations of the gravitational 
lens systems used in this paper: the PG sample and the 
JVAS/CLASS sample. 



4 CALCULATIONS 

All calculations here implement the method of PG, which 
uses the Spearman rank correlation test to generate a rel- 
ative probability for a given cosmological model. PG noted 
the fact that they always obtained a low probability with 
their sample, even when allowing for non-flat cosmological 
models (albeit in a limited area of parameter space) , galaxy 
evolution or departure from the singular isothermal sphere 
model. As PG noted, allowing for these effects increases the 
probability, since they all tend to create a negative correla- 
tion in a flat universe, but the magnitude of the effect is not 
large enough to explain the observations. Again as noted by 
PG, if the lenses are parts of clusters, then this will work in 
the opposite direction, making the observed negative corre- 
lation even more puzzling. 

Calculations were done for four samples: 

the PG sample 

the PG sample with the addition of the system 0218 + 357 
the JVAS/CLASS sample 
the union of all samples 

Note that the source redshift for 0218 + 3 57 had been pub - 
lished before the PG analysis was done ( Lawrence 1996 ). 
Since 0218 + 357 lies below and to the left of all other data 
points, it is clear that including it will weaken the puzzling 
negative correlation found by PG; this is discussed more 
quantitatively in Sect. ^. 



5 RESULTS AND DISCUSSION 

Since the PG test assigns a low probability to a k — 
universe, the question arises as to whether it can be used 
as a general cosmological test to determine the values of 
Ao and Slo- This is not the case. For all four samples I 
have calculated the Spearman rank correlation probability 
as a function of Ao and Slo in a range of parameter space 
(—8 < Ao < 2 and < fio < 10) much larger than that 
allowed even by a generous interpretation of observations. 
This was done with a resolution of 0.1 in both Ao and fio 
for both r] = 1 and rj — 0. The Spearman rank correlation 
probability is essentially constant over a wide range of pa- 
rameter space; basically, either all cosmological models are 
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probable, or all are improbable, depending on the sample 
used. 

The probability is a weak function of the cosmological 
model, with the sharpest transition occuring when crossing 
the k = line in the Ao-fio plane. For all samples except 
the PG sample, the probability is >5% in almost the entire 
parameter space;^] those cosmological models with a lower 
probability are among those ruled out by current observa- 
tions. Thus, the Spearman rank correlation probability does 
not allow one to reject any otherwise viable cosmological 
models, which shows both that there is no reason to expect 
unknown effects in the gravitational lens samples and that it 
is not very useful as a cosmological test. For the PG sample, 
the 1% contour corresponds almost exactly to the k = line, 
with higher values for a negatively curved universe. Thus, 
the PG sample is marginally compatible with a k — — 1 
cosmological model, although the probability values are low 
throughout the Ao-fio plane, with values near the maximum 
of 0.025 being attained only for small (but realistic) Qo val- 
ues and large (in absolute value) negative values of Ao- Since 
there are no known selection effects which can account for 
the differences between the PG sample and other samples, 
either the test is not very useful and/or it is pointing to un- 
known selection effects in the literature sample used by PG. 
The fact that the PG result changes dramatically (probabil- 
ity pa 10-20% in most of the Ao-fio plane) by the inclusion of 
just 1 additional data point, which could have been included 
in their analysis, argues in favour of the former possibility. 

The above discussion was for rj = 1. For rj — the sit- 
uation is qualitatively the same and quantitatively involves 
only slightly different values of probabilities derived from 
the Spearman rank correlation test. 

It is interesting to compare the probabilities from the 
Spearman rank correlation test for the PG sample using the 
actual values of z s and AO as used by PG to those obtained 
using more up-to-date data for the same lens systems. If two 
values are very near each other, rounding them off to the 
same values produces a different result for the rank correla- 
tion test than if they differ by even a small amount. Using 
more up-to-date data, an even lower probability is obtained 
for the PG sample, for r) = 1 and rj = 0, for a wide variety 
of cosmological models. 



6 CONCLUSIONS 



Park & Gott (1997) pointed out that the image separations 
in gravitational lens systems show a strong significant neg- 
ative correlation with the source redshift, while in a flat 
universe one would expect no correlation (while a negative 
correlation would be expected in a universe with negative 
curvature and a positive one in a universe of positive curva- 
ture). None of the possibilities they examined were strong 
enough to explain the effect. A possibility not examined by 
them, namely an inhomogeneous universe, produces a nega- 
tive correlation regardless of the sign of the curvature, but it 



t For the JVAS/CLASS sample, the maximal probability is 0.955 
and is realized in almost the entire k = +1 area of the parameter 
space. 



too is not strong enough to account for the effect. As a gen- 
eral test for the values of Ao and Qo the test is of no use, all 
cosmological models being assigned roughly the same prob- 
ability, but which value they are assigned depends on the 
sample used. 

The strong dependence of the result on the sample used 



seems to indicate that the result of Park & Gott (1997 



is due not to some physical cause but rather to unidenti- 
fied selection effects in the sample of gravitational lenses 
taken from the literature. The large number of JVAS and 
CLASS lenses gives us an independent comparison sample, 
thus demonstrating the need for discov ering a large numbe r 
of lenses in a well-defined sample. As Park & Gott (1997) 
point out, since many conclusions based on 'conventional' 
gravitational lensing statistics are based on essentially the 
same lenses as in their literature sample, if this sample is 
for some unknown reason atypical, then conclusions drawn 
from statistical analyses of it must be examined with care. 
It will thus be interesting to see what conclusions can be 
drawn from a statistical analysis of the JVAS/CLASS sam- 
ple after the observational tasks have been completed. (We 
expect to find more lenses, but have no qualms about using 
the present incomplete sample in this analysis since there 
is no reason to believe that a larger sample would show a 
different A8-z s relation.) 



7 NOTE 



Since this work was completed, two other responses to Park 



fe Gott (1997|) (apart fro m [Helbig (1998[ )) have appeared. 
The first (Williams 1997) is complementary to this work in 



that it assumes the effect is real and explores the a strophys- 
ical consequences while the second ( Cooray 199£ ) is more 
similar to this analysis, arriving at essentially the same con- 
clusions though using different observational data (and ex- 
ploring neither the question of usefulness as a general test 
for Ao and fio nor the effects of a locally inhomogeneous 
universe) . 
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